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The Jahn- Teller system _E (g) 61 © 62 has a particular degen- 
eracy, where the vibronic potential has an elliptical minimum. 
In the general case where the ellipse does not reduce to a cir- 
cle, the classical motion in the potential is chaotic, tending 
to trapping near one of the extrema of the ellipse. In the 
quantum problem, the motion consists of correlated tunnel- 
ing from one extremum to the opposite, leading to an average 
angular momentum reminiscent of that of the better known 
i? eg) e dynamic Jahn- Teller system. 



In the well-known E ® e Jahn- Teller (JT) effect, a 
molecule has a two-fold electronic degeneracy coupled to 
a doubly degenerate vibrational mode. This leads to a 
'conical intersection' in the vibronic potential which has a 
degenerate, circular minimum ('Mexican hat potential'), 
although higher-order vibronic coupling can break the 
ring up into three degenerate minima along the trough 
('tricorn potential') Quantum mechanically, the cou- 
pled electron-molecular vibrational (vibronic) wave func- 
tion can tunnel between the three minima, leading to a 
ground state with a net angular momentum Q . Remark- 
ably, this 'orbital' angular momentum is quantized in 
half-integer multiples of h, indicating the strong coupling 
of electronic and molecular motions. This quantization is 
a signature of the Berry phase [^IJ] of tt associated with 
the dynamic Jahn- Teller effect; the tt Berry phase has 
been experimentally verified in triangular Naa molecules 

. Points of conical intersection lead to chaotic behavior 
in the vibrational spectra, manifested quantum mechani- 
cally by anomalous level statistics However, the high 
symmetry of the E ® e problem precludes chaos , so 
multimode interactions must be included, and the chaos 
generally appears at high energies (above the conical in- 
tersection) where many vibrational modes are excited. 

Here, it is shown that a simple modification of the sym- 
metry preserves the anomalous Berry phase, yet leads 
to chaotic behavior at much lower energies, without the 
need of additional mode coupling. This case is the square 
X4 molecule with square planar symmetry, Di^h, corre- 
sponding to an E®bi®b2 Jahn- Teller problem [Q, Fig.|^. 
The high symmetry allows two JT modes, with indepen- 
dent frequencies oji, i = 1,2, and electron- vibration cou- 
plings Vi. In the special case cui = uj2, Vi — V2, the 
problem reduces exactly to that of the E (g) e molecule. 
However, there is an intermediate case, which seems not 
to have been explored till now. When Vi/uji = V2/UJ2, 
the two modes have the same JT stabilization energy. 



-^JT ~ f^-^^h ^'^d hence the vibronic potential has 
an elliptic minimum, which is not circular unless uji = uj2- 
Given the elliptic minimum, the possibility of a periodic 
orbit arises. However, angular momentum is not con- 
served. In the present paper I analyze the resulting mo- 
tion. 
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FIG. 1. Bi (a) and B2 (b) distortions of a square X4 
molecule. 

The phonon modes of amplitude Qi are defined as 
follows. The four atomic positions. Fig. |l|a, can be writ- 
ten as 

n = no - Qiy + Q2X, 
r2 = r2o - Qix - Q2y, 
rs = rso + Qiy - Q2X, 
ri = ^40 + Qix + Q2y, (1) 

where the rio's are the positions of the X atoms in the 
undistorted square. The vibronic interaction Hamilto- 
nian is 



Hvib — ViQiTx + V2Q2TJ/ 
1 \ . .r^ fO I 



-II ■ ^202, ^ Q 



(2) 



Here the electronic operators are represented by the pseu- 
dospin T!i's and other factors are included in the electron- 
phonon coupling Vi. To the vibronic Hamiltonian must 
be added an electronic term H^i and a phononic part 
Hph, with 

Hph = ^ (Pi + Pi + M^ojIQI + M^LolQl) , (3) 

with uJi the bare phonon frequencies. A spin-orbit cou- 
pling can be included ||] 



Hso — AL • S . 



(4) 
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For a static JT effect, the momenta Pi can be ne- 
glected, and the Qi are chosen to minimize the energy, 
Eqs. |,| The solution can be written in terms of the JT 



energy E)^ - 



V^/{2Mujf ). For E'^j^ ^ E%, the low- 



the 



est energy state consists of a distortion of the mode with 
larger JT energy only. For instance, if E 
solution is Qi = 0, (52 = V2/{MujI), E 
Special cases arise when 



JT ^ ^JT) 



-E 



(2) 
JT- 



'JT ~ 



p(2) = 
'JT 



(5) 



If additionally uji = lu2, the problem reduces to the 
E ® e problem, and in Eq. ^ 7 = 9/2: the electronic 
wave function is double valued: when 6 changes by 27r, 7 
has only changed by tt (the wave functions have changed 
sign) . This sign change is the signature of a Berry phase 
[^,|| , and causes the vibronic orbital angular momentum 
to take on half-integer values |Q. This can be seen as 
follows . The z-component of orbital angular momen- 
tum is Lz = (Q1-P2 — Q2Pi)/^, and the operator which 
commutes with the vibronic hamiltonian, Eq. |^ is 



Eliminating the electrons produces the vibronic potential 
surfaces 



M 



A2 
4 ' 



(6) 



When Eq. |^ is satisfied, the lower vibronic surface has a 
minimum which is degenerate along a trough, similar to 
the Mexican hat: 



Qi ^ Qo 



COS 9 



sin 61 

UJ2 



(7) 



jz — Lz 



-Tz. 



(10) 



Since Lz is quantized in integers, j2 has half-integral 
quanta. Note from Fig. I that even when toi ^ lu2, 9 must 
change by An to produce a 27r change in 7, suggesting a 
similar Berry phase. This can be directly demonstrated. 
The Berry phase is M 



= -Psl ^ 

lo 1 + sinh 



= -27rs, 



(11) 



with Qj] = y/2EjT - X'^/SEjT, and 9 arbitrary. Near 
the trough, the lower potential surface can be expanded: 



E- = [g-q) 



(8) 



with q = (gi,92), Qi ^ Qi ~ Qi, a ^ I - \^ /l&E^j, and 
g — (wi cos 9, LU2 sin 9) - that is, there is a restoring force 
only 'perpendicular' to the trough. Defining (3^^ = LU2/UJ1, 
the electronic eigenvectors are 



"0+ = COS 7-01 + sin7-!/'2 

1p- = — SUlJTpl + COS 7-02, 



(9) 



where s is half an odd integer, introduced to make the 
total wave function single valued. Thus the Berry phase 
is TT, modulo 27r, for any anisotropy. 

While this is a standard JT problem, I have not found 
any detailed analysis of the limit Eq. ||. As a first step, I 
perform a canonical transformation 



He 



-iS 



H + i[S, H] 



with 



S = -{4PiT. + -^P2Ty) 



(12) 



(13) 



where tan 7 = + Ssin^ 9 — cos9)/{PujSm6) and 6 = 
— 1, Fig.|[ By convention, UI2 is assumed to be the 
higher frequency {Pui > 1)- 




FIG. 2. Electronic phase 7 vs phononic phase 8, for 5 — 
0.1 (solid line), 1 (dotdashed line), 10 (dashed line), and 100 
(dotted Una). 



The canonical transformation can be performed exactly 
I p^ , but for present purposes only the first order result 
is needed. S, Eq. |l^, was chosen to exactly cancel the 
term linear in Q. It yields a correction 



H2 — i[S, Hyib] 



V1V2 



AT,, 



(14) 



with Tq the identity matrix. 



A = UJIP1Q2 - UJIP2Q1 = -^\Lz + ujI{PiQ2 + P2Q1), 

(15) 



and Wj. 



jI ± Lol)/2. Thus, when 



0, 



H'2 is 



proportional to LyTz-, and the angular momentum jz — 
Lz -f T2/2 is conserved (Eq. |l^). For the present case 
til- ^ and jz is not constant. 

Given the presence of a circular trough in the potential, 
circulating orbits should be possible: could it be that 
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there is a nonvanishing average < jz >7^ even though 
jz is not constant? This possibihty can be explored in 
the related classical Hamiltonian (particle in a non-linear 
potential well) by numerically integrating the equations 
of motion 



2E 



JT 



-aujfQiil 



90 



A2/4 



where the last form utilizes the quadratic approxima- 
tion, Eq. |8[ dots indicate time derivatives, and = 
2Ejtol. The integral is evaluated using a Runge-Kutta 

routine with initial conditions Q(0) = (go /wi , 0) , (3(0) = 
(0, Pqo/LU2)- In the remaining analysis, I take A = 0. 



(16) 




These values may most easily be found by plotting 4>{T) 
vs (3 for some long time T. Typical examples are il- 
lustrated in Fig. H, while the time series are shown in 
Fig. ^. Poincare maps (plots of Qi vs Qi when Q2 = 0), 
Figure |^, confirm the chaotic nature. (Note that the 
curve (7,0.05394) is almost periodic - see particularly 
Qi(t), Fig. ^ - but the Poincare map is clearly chaotic, 
Fig. |d.) While the E®e limit, /3„ = 1, IS quasiperi- 
odic (the Poincare map is a smooth closed curve), for 
/3ij 7^ 1 even the special values are weakly chaotic, with 
the Poincare maps. Fig. |^a, having a finite spread away 
from smooth curves. The similarity of these special tra- 
jectories to scars in e.g., Sinai stadia jl^] should be noted. 




FIG. 3. Winding angle (ft calculated from Eq. |l^, for sev- 
eral values of {(5^,(5). Inset = time series, Qi,Q2{t) for 
/3„ =u;2/oJi = 2, l3 = 0.5. 

Given Qi{t), Q2{t), a winding angle (j) is defined such 
that 



Q1Q2 — Q2Q1 



(17) 



If one applies this procedure to the E ® e problem 
{lj2 — uji), the results are quite simple (long dashed line 
in Fig. 1^): (f) increases linearly with time, although the 
frequency is not constant, but varies approximately log- 
arithmically with the velocity parameter /3. By contrast, 
when UJ2 7^ ldi Figure |^ shows that (j) is generically a ran- 
dom function of time, with no linearly increasing part 
indicative of a non-zero < jz >■ The various data sets 
are characterized by the two parameters (102 /uJi, (3). (The 
figure utilizes the exact form of Eq. the approximate 
form yields equivalent results.) The figure also clearly 
suggests that the motion is chaotic. This is further indi- 
cated by the direct time series, inset of Fig. |[ 

On the other hand, there are certain special values of 
the initial conditions for which the motion is approxi- 
mately periodic, and does increase linearly with time. 



FIG. 4. Time series, Q2{t) vs Qi{t) (or vs t, in (d)) for 
several choices of /3^,/3): (a) = (2,1), (b) = (3,0.75), (c,d) 
= (7,0.05394). In frames (a-c), the ellipses are equipotential 
contours, with the beaded contour representing the potential 
minimum. 
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FIG. 5. Poincare maps for (/3^,/3) = (2,1) (a), (7,0.05394) 
(b), (2,0.5) (c), (7,1) (d). In (a), one attractor is shown blown 
up. 

How is this chaotic behavior manifested in the quan- 
tum hmit? To explore this, it is convenient to first 
rescale the variables, so that the potential has circu- 
lar symmetry, and the anisotropy appears in the ionic 
mass, nii = M{ujo/uJi)'^ , with ujq = {ujf + a;|)/2, and 
then reduce the problem to one dimension by assum- 
ing that the motion is confined to the bottom of the 
trough and only varies. The Hamiltonian becomes 
H = —'h?h/{2m+p'^), where po is the equillibrium trough 
radius, = {m^^ ± mj^^)/2 and 



h = d'i+ 6![cos20(- — dj,) + 3 sin 2090] — A4 cos' 



(18) 



with a — m+/m^ — (/J^ — ^)/{0t + 1) higher or- 
der vibronic effects are incorporated in the term pro- 
portional to A4. Schroedinger's equation can be inte- 
grated numerically, letting tp{(/),t) = ip{je,nS) = 'ip", 
with = (V'^+i - and pi 



1 



or finally (1 



Tl+l 



1 + 

= (1 



(19) 
with 7 — 
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probability spreading over two adjacent wells, as shown 
in Fig. at times 4, 7, and 11. The tunneling is co- 
herent, so there is a net circulation. Additional informa- 
tion can be found by analyzing the Husimi density [T^ ] 
Ph{p, q)^\ <P, q\^p > with 



exp[ 



2h 



■P , 



(20) 



which describes the approximate smearing of ■(/; in g and 
p as a function of time. Typical results are shown in Fig. 
0, for squeezing parameter s = 1. 

Thus, the quantum system shows a 'memory' of the 
classical chaos, in that the wave function shows similar 
trapping near the points Q2 = 0. However, whereas the 
wavefunctions appear to vary stochastically from cycle to 
cycle. Fig. the average of the wave function progresses 
smoothly. Fig. |6^. The main difference is that classically, 
the wave function can be reflected from a trapping region, 
reversing its direction of motion, while the quantum wave 
function always moves in the same direction, similar to 
the classical problem with special initial conditions. It 
seems plausible to interpret the special choice of initial 
conditions as analogous to a Bohr-Sommerfeld quantiza- 
tion condition in the quantum problem. 
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FIG. 6. Quantum time evolution, showing (a) position of 
wave function peak as a function of time, and (b) actual 
distribution of \ip\^ at several equally-spaced time intervals 
(a = 0.3, Aa = 0). 

Equation ^ was integrated numerically, assuming an 
initial gaussian distribution. Figure ^jb shows \'4'{<i),t)\'^ 
for a variety of times t. The data can be better under- 
stood from Fig. |^a, which plots (t>max vs t, where (jymax 
is that value of for which l-^p has its maximum value. 
The wave function remains trapped in one of the effec- 
tive potential wells, then quickly hops to the next one in 
a relatively short time. This hopping takes place by the 



FIG. 7. Contour plot of Husimi distribution pH of data 
similar to that of Fig. ^at several time intervals. An interwell 
hopping event occurs between times 90 and 120. 

As shown in Fig. |^a, the position of the wave function 
peak has a step-like component superposed on an aver- 
age shift with time. This average shift is independent 
of the mass anisotropy |T5| , hence corresponding to the 
same quantized angular momentum as in the isotropic 
case. This is consistent with the Berry phase remaining 
TT, Eq. but somewhat surprising in light of the classi- 
cal chaos. Blumel has suggested that this might be 
a manifestation of quantum localization in angular mo- 
mentum space u%, while the classical problem leads to 
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angular momentum space diffusion. This possibility will 
be explored in future work. 

In addition to its molecular interest, the present re- 
sults may have condensed matter applications. Berryonic 
matter has been postulated to explain anomalous 
properties of Buckyballs and other dynamic JT systems, 
but based on unit cells of triatomic molecules. Poten- 
tial applications are greatly expanded for bases of square 
molecules [ p^ . 

I thank J. Jose, F.S. Ham, R. Englman, R. Blumel, 
and B. Barbiellini for stimulating conversations. 
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